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ABSTRACT 


Five coronaviruses (CoVs) were detected in diarrheal feces from four zoo ruminant 
species: one wisent (Bison bonasus), two Himalayan tahr (Hemitragus jemlahicus), one 
sitatunga (Tragelaphus spekii), and one nyala (Tragelaphus angasii). We sequenced and 
analyzed the spike (S) and hemagglutinin/esterase (HE) genes of these viruses and 
compared the nucleotide (nt) and deduced amino acid (aa) sequences with those of other 
bovine CoV (BcoV) strains. Comparison of the entire deduced aa sequences of the S and HE 
glycoproteins revealed no specific differences that would account for discrimination 
between bovine-like CoV strains from zoo ruminants and BcoVs strains. In addition, the 
99.9% aa identity among the five CoV strains revealed that the ruminants were infected by 
the same strain. Phylogenetically, bovine-like CoVs belong to group 2a CoVs, which are 
related most closely to the BcoV strains recently isolated in Korea. These data suggest that 


cattle are potential reservoirs for CoVs that are capable of infecting zoo ruminants. 


1. Introduction 


Coronaviruses (CoVs) are a member of the order 
Nidovirales and family Coronaviridae, each member of 
which contains a single-stranded, non-segmented RNA 
positive-sense genome that is 26-31 kb in size. CoVs are 
classified into three different genetic and antigenic genera, 
each of which is divided into several subgenera. Bovine- 
like CoVs and bovine CoVs (BCoVs), which belong to 
Betacoronavirus 1, contain five major structural proteins: 
nucleocapsid (N), transmembrane (M), spike (S), small 
membrane (E), and hemagglutinin/esterase (HE) (Lai and 
Cavanagh, 1997). Betacoronavirus 1 also include the human 
CoVs-OC43 (HCoV-OC43), porcine hemagglutinating 
encephalomyelitis virus (PHEV), equine coronavirus 


* Corresponding author at: Viral Disease diagnostic Lab, Animal Disease 
Diagnostic Center Division, National Veterinary Research and Quarantine 
Service, 480 Anyang-6-Dong, Anyang, 430-824, Republic of Korea. 

Tel.: +82 31 467 4344; fax: +82 31 467 1868. 
E-mail addresses: jku0623@korea.kr, jku0622@nvrqs.go.kr (J.-K. Oem). 


0378-1135/$ - see front matter © 2010 Elsevier B.V. All rights reserved. 


doi:10.1016/j.vetmic.2010.08.035 


© 2010 Elsevier B.V. All rights reserved. 


(ECoV), canine respiratory coronavirus (CRCoV) and buba- 
line coronavirus (BuCoV) (Decaro et al., 2008; Erles et al., 
2007; Guy et al., 2000; Haring and Perlman, 2001; Vijgen 
et al., 2006). 

Bovine-like CoVs were first reported in ruminant species 
in the late 1970s (Tzipori et al., 1978). Later, bovine-like 
CoVs were also identified in wild animals (sitatunga, 
waterbuck, musk oxen, llamas, alpacas, white-tailed deer, 
giraffe, Sambar deer, and sable antelope) (Cebra et al., 2003; 
Chasey et al., 1984; Hasoksuz et al., 2007; Tsumemitsu et al., 
1995). BcoVs cause newborn calf diarrhea (CD), winter 
dysentery (WD) in adult cattle, and respiratory tract illness 
in calves and feedlot cattle (Lathrop et al., 2000; Saif et al., 
1991; Snodgrass et al., 1986). BcoVs also infect the 
respiratory tract of cattle and cause disease therein (Cho 
et al., 2001; Hasoksuz et al., 1999). 

The genome of Betacoronavirus 1 encodes a short, spike- 
like protein known as HE. BcoV uses N-acetyl-9-O-acet- 
ylneuraminic acid (sialic acid) as a receptor to initiate 
infection (Schultze and Herrler, 1992). Although S protein 
has an affinity for 9-O-acetylated sialic acid, HE protein has 
been identified as the principal BcoV sialic acid-binding 
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protein (Schultze et al., 1991a,b). Therefore, HE may be 
involved in virus entry and/or virus release from infected 
cells (Schultze et al., 1991a). The S protein, which is 
composed of two subunits, S1 and S2, has several important 
functions during the virus—host interaction. The $1 subunit 
is associated with binding to host cell receptors, whereas the 
S2 subunit is the viral M protein that mediates the fusion of 
viral and cellular membranes. Thus, S protein is critical for 
virus entry and pathogenesis (Schultze et al., 1991b). The 
variation in host range and tissue tropism among CoV is 
largely attributable to variations in the S protein (Gallagher 
and Buchmeier, 2001). 

Several Korean BcoV strains have recently been charac- 
terized at the molecular level (Jeong et al., 2005a,b; Park 
et al., 2006). However, such analyses of bovine-like CoV 
Strains circulating in Korea have not yet been conducted. In 
this study, we report sequence and phylogenetic analyses of 
the S and HE protein genes of five bovine-like CoVs isolated 
from four species of zoo ruminants. 


2. Materials and methods 
2.1. Clinical signs 


The outbreak began with severe diarrhea in a wisent 
(Bison bonasus) and then spread to other ruminants in the 
National Zoo during the 2010 winter dysentery (WD) 
outbreak. Most zoo ruminants presented, apart from the 
bloody diarrhea, with weakness, depression, anorexia, and 
dehydration. Sitatunga (Tragelaphus spekii) showed both 
severe diarrhea and respiratory symptoms. Out of 16 
affected animals, a total of 12 fecal samples were collected 
from four wisents (B. bonasus), three Himalayan tahr 
(Hemitragus jemlahicus), three sitatunga (T. spekii), and two 
nyala (Tragelaphus angasii) using rectal swabs. Clinical 
signs were observed in adult ruminants (more than 24 
months of age). No cattle farms were present in the vicinity 
of the outbreak and no cattle were housed in the Zoo. All 
fecal samples were examined for common bacterial enteric 
pathogens (Salmonella spp., Clostridium spp., Campylobac- 
ter spp., Mycobacterium spp.) and parasites (Coccidium spp., 
Crytosporidium spp.). All fecal samples were also analyzed 
for the presence of BcoV, bovine rotavirus and pestiviruses 
by RT-PCR and electron microscopy (EM). No enteric 
pathogens excepting BcoV were detected and isolated in all 
fecal samples. 


2.2. RNA extraction and RT-PCR 


Total RNA was extracted from clarified fecal samples 
using an RNeasy minikit (Qiagen, Valencia, CA) according 
to the manufacturer’s directions. One-step RT-PCR reac- 
tions were performed as previously described (Jeong et al., 
2005b; Park et al., 2006). 


2.3. Sequence determination 


PCR products were purified using an agarose gel DNA 
extraction kit (INtRON Biotechnology Inc., Daejeon, Korea) 
and subcloned into the pGEM-T vector (Promega Corp., 
Madison, WI), according to the manufacturer’s directions. 


Automated nucleotide sequencing of the 5’-UTR gene, 
inserted into the vector, was performed using a genetic 
analyzer and cycle sequencing kit (Applied Biosystems, 
Foster City, CA). All nucleotide positions were confirmed by 
three or more independent sequencing reactions in both 
directions. 


2.4. Sequence comparison and phylogenetic analysis 


Nucleotide sequence alignments were generated using 
the BioEdit computer program (Ibis Biosciences Inc., 
Carlsbad, CA). For phylogenetic analyses, complete 
sequences of the S and HE genes from four zoo ruminants 
were aligned with reference sequences retrieved from 
GenBank (Table 1). A phylogenetic analysis based on the S 
and HE sequences was also conducted using BioEdit with 
the Molecular Evolutionary Genetics Analysis (MEGA) 3.1 
program, with bootstrap values based on 3000 replicates. 
The cutoff point for bootstrap replication was >70%. 


3. Results and discussion 


Of the 12 fecal samples tested, five (one wisent, two 
Himalayan tahr, one sitatunga, and one nyala) were 
positive for CoV. Five fecal samples were used for virus 
isolation and sequencing. All of the S genes sequenced 
contained an open reading frame (ORF) of 4092 nucleo- 
tides, encoding a predicted protein of 1363 amino acid 
(aa) residues and molecular weight of approximately 
150 kDa. Sequence analysis of the five identified strains 
Showed a 99.9% amino acid identity and a genetic 
relatedness of 99.4-99.5% to two other Korean isolates, 
BcoV-0501 and BcoV-0502. Thus, the most recent 
isolates, BcoV-0501 and BcoV-0502 (identified in 
2005), had the highest nucleotide identity with the five 
bovine-like CoVs, whereas they showed only 97.8-97.9% 
identities with the BC94/1994 (DQ389641) and Sun5/ 
1994 (DQ389657) strains. In total, 86 nucleotide sub- 
stitutions were found compared to the Mebus strain, but 
no insertions or deletions were detected in any sequence. 
Also, 34 substitutions were found in the complete aa 
sequence of the viral glycoprotein S. 

A phylogenetic analysis using reference strains from all 
known groups of CoVs demonstrated that the bovine-like 
CoVs isolated in this study belonged to Betacoronavirus 1 
clustering with Korean BcoVs isolates collected after 2004 
(Fig. 1a). The bovine-like CoVs isolated were closely related 
to two Korean strains (BcoV-0501 and BcoV-0502) identified 
in 2005, suggesting that this virus may have been transmit- 
ted from cattle to zoo ruminants. Korean BcoVs strains 
collected before 2004 were divided into several clusters. A 
BcoV-white-tailed deer strain isolated in 1994 (Tsumemitsu 
et al., 1995) formed a cluster together with Korean BcoV- 
KWD1-10 isolates, with the exception of KWD 7. 

All of the HE genes contained an ORF of 1272 nucleotides, 
encoding a predicted protein of 424aa residues and 
molecular weight of approximately 47.6 kDa. In total, 27 
nucleotide substitutions were found compared to the 
Mebus strain, but no insertions or deletions were detected 
in any sequence. Also, five substitutions were found in the 
complete aa sequence of the S glycoprotein. 
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Table 1 

GenBank accession numbers of Betacoronavirus 1 reference strains used in the phylogenetic analysis. 
Strain Year Host 
Wisent/2010 2010 Wisent 
Sitatunga/2010 2010 Sitatunga 
Nyala/2010 2010 Nyala 
Himalayan Tahr 1/2010 2010 Himalayan tahr 
Himalayan Tahr 2/2010 2010 Himalayan tahr 
0501/2005 2005 Cattle 
0502/2005 2005 Cattle 
BC94/1994 1994 Cattle 
SUN5/1994 1994 Cattle 
A3/1994 1994 Cattle 
KCD1/2004 2004 Calf 
KCD2/2004 2004 Calf 
KCD3/2004 2004 Calf 
KCD4/2004 2004 Calf 
KCD5/2004 2004 Calf 
KCD6/2004 2004 Calf 
KCD7/2004 2004 Calf 
KCD8/2004 2004 Calf 
KCD9/2004 2004 Calf 
KCD10/2004 2004 Calf 
KWD1/2002 2002 Cattle 
KWD2/2002 2002 Cattle 
KWD3/2002 2002 Cattle 
KWD4/2002 2002 Cattle 
KWD5/2002 2002 Cattle 
KWD6/2002 2002 Cattle 
KWD7/2002 2002 Cattle 
KWD8/2002 2002 Cattle 
KWD9/2002 2002 Cattle 
KWD10/2002 2002 Cattle 
KWD11/2002 2002 Cattle 
KWD12/2002 2002 Cattle 
KWD13/2002 2002 Cattle 
KWD14/2002 2002 Cattle 
KWD15/2002 2002 Cattle 
KWD16/2002 2002 Cattle 
KWD17/2002 2002 Cattle 
KWD18/2002 2002 Cattle 
KWD19/2002 2002 Cattle 
White-tailed deer/1994 1994 White-tailed deer 
Calf-giraffe/2003 2003 Giraffe 
Alpaca/1998 1998 Alpaca 
LUN/1998 1998 Cattle 
ENT/1998 1998 Cattle 
Sable antelope/2003 2003 Sable antelope 
Sambar deer/1994 1994 Sambar deer 
Waterbuck/1994 1994 Waterbuck 
LSU/1994 1994 Cattle 
OK/1996 1996 Cattle 
E-DB2/1984 1984 Cattle 
LY-381 1965 Cattle 
Nordern vaccine/1991 1991 Cattle 
L9/1991 1991 Cattle 
Mebus/1972 1972 Cattle 
Quebec/1972 1972 Cattle 
Bubalus/2007 2007 Water buffalo 


Clinical forms Country GenBank accession no. 

S HE 
Winter dysentery Korea HM573326 HM573321 
Winter dysentery Korea HM573329 HM573324 
Winter dysentery Korea HM573330 HM573325 
Winter dysentery Korea HM573327 HM573322 
Winter dysentery Korea HM573328 HM573323 
- Korea EU686689 EU401975 
- Korea EU401986 EU401676 
- Korea DQ389641 EU401979 
- Korea DQ389657 EU401978 
- Korea AY935640 EU401977 
Calf diarrhea Korea DQ389632 DQ389642 
Calf diarrhea Korea DQ389633 DQ389643 
Calf diarrhea Korea DQ389634 DQ389644 
Calf diarrhea Korea DQ389635 DQ389645 
Calf diarrhea Korea DQ389636 DQ389646 
Calf diarrhea Korea DQ389637 DQ389647 
Calf diarrhea Korea DQ389638 DQ389648 
Calf diarrhea Korea DQ389639 DQ389649 
Calf diarrhea Korea DQ389640 DQ389650 
Calf diarrhea Korea DQ389641 DQ389651 
Winter dysentery Korea AY935637 DQ016118 
Winter dysentery Korea AY935638 DQ016119 
Winter dysentery Korea AY935639 DQ016120 
Winter dysentery Korea AY935640 DQ016121 
Winter dysentery Korea AY935641 DQ016122 
Winter dysentery Korea AY935642 DQ016123 
Winter dysentery Korea AY935643 DQ016124 
Winter dysentery Korea AY935644 DQ016125 
Winter dysentery Korea AY935645 DQ016126 
Winter dysentery Korea AY935646 DQ016127 
Winter dysentery Korea DQ389652 DQ994162 
Winter dysentery Korea DQ389653 DQ994163 
Winter dysentery Korea DQ994164 DQ389654 
Winter dysentery Korea DQ994165 DQ389655 
Winter dysentery Korea DQ994166 DQ389656 
Winter dysentery Korea DQ994167 DQ389657 
Winter dysentery Korea DQ994168 DQ389658 
Winter dysentery Korea DQ994169 DQ389659 
Winter dysentery Korea DQ994170 DQ389660 
Sporadic diarrhea USA FJ425187 FJ425187 
Winter dysentery USA EF424623 EF424623 
Severe diarrhea USA DQ915164 DQ915164 
Bovine shipping fever USA AF391542 AF391542 
Bovine shipping fever USA AF391541 AF391541 
Diarrhea USA EF424621 EF424621 
Winter dysentery USA FJ425188 FJ425188 
Winter dysentery USA FJ425186 FJ425186 
Respiratory problem USA AF058943 AF058943 
Respiratory problem USA AF058944 AF058944 
- USA DQ811784 DQ811784 
- USA AF058942 AF058942 
- USA M64668 - 
- USA M64667 M76327 
- USA U00735 U00735 
- Canada AF220295 - 
Severe diarrhea Italy EU019216 EU019216 


Phylogenetic analyses were carried out to determine 
genetic relationships based on the nucleotide sequences of 
the HE gene (Fig. 1b). The Korean BcoV strains were divided 
into three groups. The first group consisted mainly of all of 
the Korean isolates. CoV strains from zoo ruminants were 
grouped within a cluster of other Korean CoV strains 
isolated after 2004. The CoV strains clustered with the 
BcoV-0501/2005 strain, which was identified in 2005. The 


second group included some Korean strains isolated in 
2002 and bovine-like strains from wild ruminants. The last 
group consisted mainly of BcoV vaccine strains. 

In this study, we compared the S and HE genes of 
Korean bovine-like CoVs isolated directly from feces of zoo 
ruminants to determine their relatedness to bovine-like 
CoVs from wild ruminants at other locations. However, no 
relationship between Korean bovine-like CoVs and those 
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(a) BcoV-Nyala/2010 
BeoV-Sitatunga/2010 

100 | |BcoV-Himalayan Tahr 2/2010 
BcoV-Himalayan Tahr 1/2010 
BcoV-Wisent/2010 
BcoV-0502/2005 


| 99 Bco V-0501/2005 
BcoV-KCD2/2004 


79 BcoV-KCD4/2004 
97 BcoV-KCD1/2004 
BcoV-KCD3/2004 
BcoV-KCD8/2004 
BcoV-KCD10/2004 
94 Bco V-KCD5/2004 
BcoV-KCD6/2004 
BcoV-KCD9/2004 
Bcov-KCD7/2004 
Bco V-KWD15/2002 


92 


99 Bco V-KWD7/2002 
BcoV-KWD18/2002 
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BcoV-KWD1 4/2002 


BcoV-KWD17/2002 
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Ld Beg V-KWD16/2002 


99 BcoV-Alpaca/1998 
L—_____ Beo V-LUN/1998 
BcoV-ENT/1998 
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BcoV-White-tailed deer/1994 
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BcoV-KWD2/2002 
i Bco V-KWD3/2002 

BcoV-KWD6/2002 
BcoV-KWD4/2002 
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BcoV-Sambar deer/1994 
1400 — BcoV-Waterbuck/1994 
Bco V-LSU/1994 
BcoV-OK/1996 


13 


Bco V-Bubalus/2 007 
BcoV-E-DB2/1984 
BcoV-LY-138/1965 
BcoV-Norden vaccine/1991 


99 BcoV-L9/1991 
100 87 Bco V-Mebus/1972 
Bco V-Quebec/1972 
77 BcoV-A3/1994 
BcoV-BC94/1994 


100 | Reo V-SUN5/1994 
0.002 


399 


Fig. 1. Phylogenetic tree generated by neighbor-joining analysis of genetic distance in the complete sequences of the S (a) and HE (b) genes. The nucleotide 
sequences of five bovine-like CoV strains were aligned using BioEdit and Mega 3 software. Three thousand bootstrap replicates were subjected to both 
nucleotide sequence distance and neighbor-joining methods. The consensus phylogenetic tree is shown with bootstrap values +>70% displayed above the 


tree branches. 


from wild ruminants at other locations was detected. 
Sequence analysis of the S and HE genes, however, 
indicated that the strains show a high similarity to BcoV. 
This was confirmed by phylogenetic analysis, thus 
supporting data obtained previously (Alekseev et al., 
2008). The CoV strain isolated in this study appears to 
derive from the same ancestor as the BcoV 05/01 and BcoV 
05/02 strains isolated from cattle in Korea. Also, the 


sequences of the S and HE genes of Korean CoV strains 
represented on the tree grouped according to year of 
isolation. These data support the hypothesis that genetic 
difference may be related more to the time of the 
appearance of an outbreak than to its geographic location 
(Alekseev et al., 2008; Jeong et al., 2005b). 

The 99.9% aa identity of the five bovine-like CoV strains 
suggests that the zoo ruminants were infected by the same 
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Fig. 1. 


bovine-like CoV. The high genetic similarity with bovine- 
like CoVs from the zoo ruminants and Korean BcoV strains 
was also confirmed, suggesting that cattle may be potential 
reservoirs for CoVs that infect zoo ruminants. However, 
additional isolation and characterization of bovine-like 
CoVs from other heterologous hosts are required to clarify 
the origin of bovine-like CoVs and their relationship to 
other BcoVs. 

Since the 1970s, bovine-like CoVs have been recognized 
as being distributed worldwide, with many reports 


BcoV-Sambar deer/1994 
BcoV-Waterbuck/1994 


BcoV- Bubalus/2007 


BcoV-LY-138/1965 
Bcov-Mebus/1972 
Bco V-Quebec/1972 
BcoV-L9/1991 


| BooV-BC94/1994 
88 Bco V-SUNS/1994 
4 'BcoV-A3/1994 


rf 


99 


(Continued ). 


indicating that the virus had been detected and/or isolated 
in the United States, Europe, and Asia. BcoVs have been 
detected in the feces of calves and adult cattle with 
diarrhea in Korea (Jeong et al., 2005a,b; Park et al., 2006). In 
addition, antibodies to BcoV have been demonstrated in 
Korean native goats (Capra hircus), and 1.0% of Korean 
native goats were serologically positive for BcoV (Yang 
et al., 2008). However, this to our knowledge is the only 
report of detection and characterization of bovine-like 
CoVs in Korea that has been published to date. No specific 
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genomic characteristics that could account for the host 
specificity exhibited by bovine-like CoVs were identified 
by comparative sequence analysis. 
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